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Abstract 

Objectives: According to the Korea Centers for Disease Control and Prevention, 
new active tuberculosis (TB) cases have Increased since 2001. Some key factors 
explain and characterize the transmission dynamics of Korean TB infection, such 
as a higher ratio of latent individuals and a new reporting system Implemented In 
2001 , among others. 

Methods: We propose a mathematical TB model that includes exogenous rein- 
fection to gain a better understanding of the recent trend for TB Incidence. We 
divide the simulation time window into two periods, 1970—2000 and 2001—2012, 
according to the Implementation date of a new TB detection system. 
Results: Two sets of parameters, including the transmission rate, the latent 
period, the recovery rate, and the proportion of exogenous reinfection, are 
estimated using the least-squares method and calibrated to data on the inci- 
dence of active TB. 

Conclusion: Among some key parameters In the model, the case finding effort 
turned out to be the most significant impacting component on the reduction In 
the active TB cases. 
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1. Introduction 

Tuberculosis is a bacterial disease caused by the bac- 
terium "Mycobacterium tuberculosis". Almost 10 million 
people are infected every year and two millions of them 
are killed by tuberculosis all over the world. Tuberculosis 
is a major cause of death for children, young adults and 
women of reproductive age, especially, in developing 
countries [1]. The Republic of Korea (hereafter "Korea") 



still has the highest infection rate and the highest mor- 
tality rate of TB infection among OECD countries in 2005 
[2], even though it has been successful to have a dramatic 
reduction in the TB incidence last few decades. Annual 
new TB incidence rate is 87 people (per 100,000 people) 
that much higher than the OECD average of 17.7. Also, 
mortality induced by active TB cases is 10 people (per 
100,000 people) that about five times higher than the 
average mortality rate 2.1. Recently, Korea Centers for 
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Disease Control and Prevention reported new active-TB 
cases have increased particularly since 2001 [3]. One of 
main reasons for sudden increase of active-TB incidences 
is due to the start of a new tuberculosis surveillance 
system implemented by the Korean government in 2000. 
Moreover, there are other important factors for explain- 
ing of the current TB transmission dynamics such as in- 
crease of aged population, a still higher proportion of the 
latent individuals after the Korean War, a higher pro- 
portion of the retreated individuals and the emergence of 
multidrug-resistance TB [4]. 

It is extremely challenging to have a full under- 
standing of the complex dynamics in any infectious 
disease. However, mathematical modeling, mathemat- 
ical analysis and numerical simulations have been 
providing a complementary tool to understand the dis- 
ease dynamics better [5]. There have been quite a 
number of researchers working on the mathematical 
modeling for TB [6—8]. Unlike the fast dynamics of 
influenza (a time scale is much faster than the one of 
TB), TB is a slow disease since it has a short and in- 
fectious period while having a long and variable latency 
period [9]. About a third of the world's population, or 
around 2 billion people, are the reservoir of the TB 
bacteria but most latently infected individuals do not 
become actively infectious (active TB). Around 10% of 
people infected with TB actually develop the disease 
within a five-year window while others do it after a long 
period of time (maybe not develop the disease at all for 
their entire life). On the other hand, infectious in- 
dividuals remain for relatively shorter periods of time 
partially because of the antibiotic treatments. 

This leads for us to take account of this long latent 
period as an important component in the TB progression 
model. The disease progression towards active TB can be 
accelerated with re-exposure to TB bacilli through 
repeated contacts with individuals with active TB 
resulting in exogenous reinfection which implies that 
latent individuals acquire a new infection from another 
infectious individual or endogenous infection, namely, 
reactivation of latent individuals [10—12]. In particular, 
the situation of the Republic of Korea such as a high 
latent ratio, exogenous reinfection should be a part of the 
model and it might play a key role in the Korean TB 
transmission. Recently, a Korean TB dynamical model 
has been developed including time-dependent co- 
efficients and further optimal control theory has been 
employed to observe what would be more effective 
control measures to minimize the number of infected and 
exposed individuals [13]. In this manuscript, we modify 
this existing model by incorporating the exogeneous 
reinfection term from the work [12]. Parameter estima- 
tion for our new TB model is carried out using the least- 
square method to the TB incidence from two periods, 
1970-2000 and 2001-2012. The model output and the TB 
incidence data are in a very good agreement. Based on 



the model and the parameters, sensitivity analysis is 
carried out in order to find out which parameters are key 
factors to reduce the active TB incidence. 



2. IMaterials and methods 

As studied in the earlier work, the standard SEIL 
compartmental model is taken as the base model for the 
TB transmission dynamics, which uses the Korean de- 
mographic parameters [13]. In addition, the exogeneous 
reinfection has been incorporated in the model (the 
mathematical model is given below). The new TB trans- 
mission model classifies individuals as the susceptible 
class (S ), the exposed class(ii ), (or high-risk latent, that 
is, recently exposed but not yet infectious), the active-TB 
infectious class (/ ), and the low risk latent class or 
treated infected class (L). 

^=PcS^~pficE^~{^i + k+a)E+crpcL^ 

dl I 

—r = jl-- (^pcL— ~ uL + aE 

dt N 

N=S+E+1 +L 

The birth and mortality rates per capita per unit of 
time are A and ^u, respectively. 13c is the transmission rate 
for the susceptible, exposed and low risk latent in- 
dividuals with an infectious individual per contact per 
unit of time. It is assumed that it is less likely to get 
infected for the low risk latent individuals with a 
reduction constant, 0<(r<l(5). k is the rate at which an 
individual leaves the latent class by becoming 
infectious, d is the per-capita disease-induced death rate 
and y is the per-capita treatment rate. Individuals who do 
not progress from the class E to the class / are moving to 
the low risk latent class at the rate a, (e.g., taking the TB 
medications before occuring active TB or called "case 
finding effort"). The term p(3cE{I/N ) models the exog- 
enous reinfection rates with p representing the level of 
reinfection. A value of 0<p< 1 implies that reinfection is 
less likely than a new infection. Most developed coun- 
tries have low incidence TB rate, therefore, the exoge- 
neous reinfection can be ignored. As mentioned earlier, 
this is not the case of Korea, since Korea has a higher 
proportion of exposed or low risk latent individuals like 
developing countries. Hence, it is better to take account 
of exogeneous reinfection in the TB transmission model. 
Also, we assume a constant per-capita removal rate to 
focus exclusively on the role of exogenous reinfection 
then, the basic reproductive number Rq for (1) can be 
obtained as follows: 
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Table 1. 
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Parameter 


Description 


Value(Tl) 


Value(T2) 


Reference 




IvdlC iiUlll Lj LU 1 Uy LUllLtlt^Llllg lO 


0 0 


0 0 


n 71 


k 


Rate from E to I 


0.005 


0.005 


ri2i 


d 


Per-captiadeasease-induced deatii rate 


0.0943 


0.0943 


[13] 


a(t) 


Non-progress rate from E to I 


0.06787 


0.0352 


Data fitted 


Pit) 


Level of reinfection 


0.2958 


0.2469 


Data fitted 


Jit) 


Treatment rate 


0.81 


1.5179 


Data fitted 


3c(0 


Number of new infectious 


6.4737 


18.673 


Data fitted 



R = ( ^ ( ^ \ 

° \ii + y + d)\iJL + k+a) 

This dimensionless quantity is represented by the 
product ficlijj, + y + d), that is, by the average popu- 
lation of infected people from one infectious individual 
during his or her infectious period and kl{ji+k+a), the 
rate of the population, which survives the exposed 
period and go to the infectious class successfully. 
Therefore Rq gives the number of secondary infectious 
cases produced by an infectious individual for a period 
of his or her infectious period in a population of sus- 
ceptible individuals. 



3. Results 

The parameter estimation of the system (1) to the 
incidence data of active-TB cases in Korea from 1970 to 
2011 is carried out. We first assume that iJL{t) and b{i) is 
time-dependent parameter for every year taken from 
[13]. The other parameters for a, k, d are defined as 
constant taken from [12]. The rest of parameters are 
estimated through the least-square method: y(?) is an 
estimate of the treatment rate; /?c(/) is an estimate of the 
transmission rate; a{f) is an estimate of the case finding 
rate from E to /; p(f) is an estimate of the level of 
reinfection. 



For the parameter estimation procedure, the least- 
square method is used (built-in-function in Matlab). We 
split up the time interval [1970, 201 1] into two intervals 
[1970, 2000], and [2001, 2011] due to the imple- 
mentation of a new tuberculosis detection system in 

2000. Now, we have estimations of yit), fic(t), a{t) and 
p{t) for two distinct time windows, which are listed in 
Table 1. Figure 1 displays the actual data of active-TB 
incidence and its best fit or the model outputs from 
1970 to 2011 yr. The active-TB incidence (shown as°) 
represents the number of new infectious individuals per 
each year compared with the model outputs, dashed and 
solid curves under two different sets of parameter values 
estimated in Table 1. The results demonstrate that our 
fitted TB model agrees well with the reported TB inci- 
dence data. 

Next, using the parameter values estimated through 
the least-square method, the model output is illustrated 
in Figure 2; the proportions of 5(0W(0, E{t)INit), 1(1)/ 
N{t), and L(t)IN(t). The fraction of susceptible in- 
dividuals is increased from aroud 0.36 to over 0.5 in a 
steady manner during the simulated time. The fraction 
of exposed individuals and the fraction of infectious 
individuals have been decreased significantly from 1970 
until 2000 but they have been slightly increased after 

2001. The fraction of latent individuals stays around 0.4 
throughout the entire period of time. When we take a 
closer look at the fraction of active TB infected 



dS I 

-=b{t)N~ficit)S-~li{t)S 

'^=p{t)pc{t)E!-+kE ~ + j{t) + d)I 
dt N 

^ = 7{t)I - a{t)fic{t)L^ - ^i{t)L + a{t)E 

'§^pmit)E^ + kE 
N=S+E+I +L 
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Figure 1. The reported incidence of active-TB (circle) and its best fit (dashed and solid curves), incidence are displayed. 



individuals, the number of infectious individuals is ex- 
pected to grow slowly since 2001. 

Again, using the same parameter values that obtained 
above, we present some sensitivity analysis to find out 
which parameters are ksy factors to reduce the active TB 
incidence. Specifically, we focus on the treatment 
effort(y), the case holding effort((r) and the case finding 
effort(a) in the model equations (1) as noted in [14]. 
Numerical simulations are performed using three 
different values for each of these effort parameters. The 
effort y denotes the effort of treatment of active TB 
infected individuals in / while a denotes the effort of 
preventing low risk latent individuals from going back to 
the high risk or exposed class by forcing them to finish 
their treatment. Lastly, a denotes the effort of identifying 
exposed individuals so that they can be treated at their 



early stage and move to the low risk class. Here only the 
sensitivity analysis for the case finding effort is shown 
since this parameter is the most significant impact on the 
reduction of active TB cases. 

Figure 3 presents the proportions of S{t)/N{t), E{t)l 
N(t), Iit)IN{t), and L{t)IN{t) using three different values 
of a to illustrate the impact of the case finding on the 
transient dynamics of TB after 2011 up to 2030. The 
corresponding results are plotted in red dotted curve, blue 
solid curve and black dashed curve using a = 0.0252 and 
a = 0.0352 (the baseline value fi-om the estimation as 
given in Table 1) and a = 0.0752, respectively. As shown 
in Figure 3, the proportion of active-TB incidences /(/)/ 
N{t) will be as twice as the one in 2011 year when 
a = 0.0352 (or even more so with a = 0.0252). When we 
put more effort on this case finding rate as a = 0.0752, the 




1970 1980 1990 2000 2010 1970 1980 1990 2000 2010 




qI . . . J o' ' ■ ' ^ 
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Time (yr) Time (yr) 

Figure 2. The proportions of S(t)/N(t), E(t)/N(t), I(t)/N(t) and L(t)/N(t) are plotted. 



44 



S. Kim, et al 



0.5 



0.4 



0.3 



2015 2020 2025 2030 



0.4 
0.3 
^0.2 
0.1 

0 





a=0.0252 




a=0.0352 




a=0.0752 











2015 2020 2025 2030 




2015 2020 2025 2030 
Time (yr) 



2015 2020 2025 2030 
Time (yr) 



Figure 3. The fractions of state variables, S(t)/N(t), E(t)/N(t), I(t)/N(t) and L(t)/N(t), are plotted. 



proportion of active-TB incidences I{t)IN{t) will be 
almost half around 2030 year. The results of sensitivity 
analysis for the other two cases, the treatment efforts and 
the case holding effort, are similar as the ones of the case 
finding effortof course, not as significant as the case 
finding case but some impact to a certain extent. 

This is consistent with the results for the basic 
reproductive number given earlier. Because the case 
holding control has no influence on basic reproductive 
number, the treatment effort and the case finding effort 
have influences on the basic reproductive number in a not 
simple way. More rigorous mathematical analysis will be 
needed for clarifying the impact of these efforts on the 
dynamics in terms of the basic reproduction number. 

4. Discussion 

Mathematical modeling in infectious diseases has 
played an important role to give us some insights on 
which counter measures would be more effective for 
reducing the negative impact of active-TB incidence. In 
this manuscript, we have modified a mathematical 
model for the TB dynamics in order to look into the role 
of exogeneous reinfection and moreover, in order to 
understand the recent increase in active-TB incidences 
in Korea. The importance of exogenous reinfections is 
still under debate. There have been some evidences, 
which support endogenous tuberculosis may be the 
main cause of tuberculosis while other using the same 
epidemiological data, has also shown the importance of 
exogenous reinfection in the USA. In the Netherlands, 
the contribution of endogenous and exogenous reinfec- 
tion to the total disease load was studied [12 and ref- 
erences are therein]. We would like to take more 
rigorous study to investigate these contributions in the 
Korean TB transmission in the near future. The 



proposed model output for the active-TB incidence 
shows a good agreement with the observed data. Among 
some key parameters in the model, the case finding 
effort turned out to be the most significant impacting 
component on the reduction in the active TB cases. 
However, our results show that concentrating on treat- 
ment alone or case finding alone will not dramatic affect 
the reduction in the active TB incidences. Therefore, 
taking two or more of the key parameters at the same 
time will go a long way in reducing the burden of the 
active TB. More systematic sensitivity analysis and 
extensive mathematical analysis should be and will be 
carried out in our future research. 
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